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ABSTRACT 

Aims. Giant radio halos are diffuse, Mpc-scale, synchrotron sources located in the central regions of galaxy clusters and provide the 
most relevant example of cluster non-thermal activity. Radio and X-ray surveys allow to investigate the statistics of radio halos and 
may contribute to constrain the origin of these sources and their evolution. 

Methods. We investigate the distribution of clusters in the plane X-ray (thermal, Lx) vs synchrotron (non-thermal, f 1 4) luminosity, 
where clusters hosting giant radio halos trace the Pia-Lx correlation and clusters without radio halos populate a region that is well 
separated from that spanned by the above correlation. The connection between radio halos and cluster mergers suggests that the clus- 
ter Mpc-scale synchrotron emission is amplified during these mergers and then suppressed when clusters become more dynamically 
relaxed. 

Results. In this context, by analysing the distribution in the Pia~Lx plane of galaxy clusters from X-ray selected samples with ade- 
quate radio follow up, we constrain the typical time-scale of evolution of diffuse radio emission in clusters and discuss the implications 
for the origin of radio halos. 

Conclusions. We conclude that cluster synchrotron emission is suppressed (and amplified) in a time-scale significantly smaller than 
1 Gyr We show that this constraint appears difficult to reconcile with the hypothesis that the halo's radio power is suppressed due 
to dissipation of magnetic field in galaxy clusters. On the other hand, in agreement with models where turbulent acceleration plays 
a role, present constraints suggest that relativistic electrons are accelerated in Mpc-scale regions, in connection with cluster mergers 
and for a time-interval of about 1 Gyr, and then they cool in a relatively small time-scale, when the hosting cluster becomes more 
dynamically relaxed. 

Key words, particle acceleration - radiation mechanisms: non-thermal - galaxies: clusters: general - radio continuum: general - 
X-rays: general 



1. Introduction 

Radio observations of galaxy clusters unveil the presence of rel- 
ativistic particles and magnetic fields in the intracluster medium 
(ICM) through the detection of diffuse Mpc-scale synchrotron 
emission in the form of radio halos and radio relics (e.g., Ferrari 
et al. 2008 for recent review). 

Giant radio halos provide the most spectacular evidence of non- 
thermal phenomena in the ICM. They are giant diffuse radio 
sources located at the centre of galaxy clusters and extending 
similarly to the hot ICM; remarkably they are always found in 
clusters with evidence for ongoing mergers (e.g. Buote 2001; 
Govoni et al. 2004; Venturi et al. 2008). These halos prove that 
mechanisms of in situ particle acceleration or injection are ac- 
tive in the ICM since the diffusion time necessary to the radio 
emitting electrons to cover Mpc scales is much longer than their 
radiative lifetime (e.g., Jaffe 1977). 

Correlations between the radio power at 1 .4 GHz of giant radio 
halos (Pi ,4) and their physical size, and between Pi 4 and the 
X-ray luminosity {Lx) and temperature of the hosting clusters 
have been found and discussed in the literature (e.g. Liang et al. 
2000; Bacchi et al. 2003; Cassano et al. 2006,07; Brunetti et al. 
2007; Rudnick et al. 2009). These correlation suggest that grav- 
ity provides the reservoir of energy to generate the non-thermal 
components responsible for the emission from the ICM. 



Mergers drive shocks and turbulence in the ICM that may 
lead to the amplification of the magnetic fields (e.g., Dolag et 
al. 2002; Subramanian et al. 2006; Ryu et al. 2008) and to the 
acceleration of high energy particles (e.g., EnBIin et al. 1998; 
Roettiger et al. 1999; Sarazin 1999; Blasi 2001; Brunetti et al. 
2001, 2004; Petrosian 2001; Fujita et al. 2003; Ryu et al. 2003; 
Hoeft & Briiggen 2007; Brunetti & Lazarian 2007). More specif- 
ically, extended and fairly regular diffuse radio emission may be 
produced by secondary electrons injected during proton-proton 
collisions, since relativistic protons can diffuse on large scales 
(hadronic or seconday models; e.g., Dennison 1980; Blasi & 
Colafrancesco 1999), or by relativistic electrons re-accelerated 
in situ by MHD turbulence generated in the ICM during cluster- 
cluster mergers (re-acceleration models; e.g., Brunetti et al. 
2001; Petrosian 2001). Observations provide support to the idea 
that turbulence may play a role in the particle acceleration pro- 
cess (e.g., Brunetti 2008; FeiTari et al. 2008; Cassano 2009 for 
recent reviews). Low frequency radio observations (e.g., with 
LOFAR, LWA) and high energy observations (with FERMI) are 
expected to set crucial constraints. 

In this paper we study the distribution of X-ray selected 
galaxy clusters in the f 1 4-Lx plane providing an extension of 
a previous work (Brunetti et al. 2007). More specifically, we dis- 
cuss constraints on the relevant time-scales of the evolution of 
magnetic fields and emitting particles in the ICM and their con- 
sequence on the origin of radio halos. 
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Fig. 1. Left panel: distribution of GMRT galaxy clusters (blue) and of other radio-halo clusters from the literature (filled black 
symbols) in the Pi 4-0.1-2.4 keV luminosity plane (Table 1). Empty circles mark giant radio halos from the GMRT sample, empty 
triangles mark the two mini-halos in cool-core clusters from the GMRT sample, the cross marks the position of RXJ1314, and 
arrows mark upper limits for GMRT clusters with no evidence of Mpc-scale radio emission. The solid line gives the best fit to the 
distribution of giant radio halos (BCES Bisector, Table 2). Right panel: distribution of giant radio halos (GMRT + literature. Table 
1) in the Pi 4-bolometric X-ray luminosity plane. The solid line gives the best fit to the distribution of giant radio halos (BCES 
Bisector, Table 2). 



In Sect. 2 we discuss the distribution of X-ray luminous 
galaxy clusters in the P\n-Lx plane and the P\n-Lx correla- 
tion traced by radio halos. In Sect. 3 we discuss the connection 
between mergers and radio halos and their evolution driven by 
these mergers. In Sects. 4 and 5 we constrain the evolution time- 
scale of radio halos and compare our results with model expec- 
tations, respectively. In Sect. 6 we give our Conclusions. 

Ho-7Q km s ' Mpc Q,„=0.3, Q.a-Q.7 are adopted throughout 
the paper 



2. The Radio - Lx correlation & cluster bi-modality 

The "GMRT Radio Halo Survey" (Venturi et al. 2007, 2008) 
has provided a significant step to a statistically solid exploration 
of the properties of radio halos through a large observational 
project carried out with the Giant Metrewave Radio Telescope 
(GMRT, Pune-India) at 610 MHz. This pointed-radio survey 
completed the radio follow up of a complete sample of 50 X- 
ray luminous {Lx > 5 • 10^"* erg/s) galaxy clusters in the red- 
shift range 0.2 - 0.4 (taken from the REFLEX, Boehringer et 
al. 2004 and the extended BCS, Ebeling et al. 1998, 2000 cata- 
logues) through high sensitivity observations of 34 clusters with 
no radio information. Large scale synchrotron emission at level 
of presently known radio halos was found only in ~ 30% of the 
selected (X-ray luminous and massive) clusters (Brunetti et al. 
2007; Venturi et al. 2008), with evidence that this fraction de- 
pends on cluster X-ray luminosity (Cassano et al. 2008). 

Figure 1 shows the distribution of GMRT galaxy clusters 
(blue) in the Pi^-Lx plane (Left : 0.1-2.4keV luminosity; Right 
: bolometric luminosity), together with that of clusters hosting 
giant radio halos (from the literature. Table 1): giant radio halos 
trace the P1.4 - Lx correlation, clusters with no large scale radio 
emission populate the region of the radio upper limits that is well 
separated from that spanned by radio halos. 



The distribution of giant halos across the correlation is signifi- 
cantly broader than the typical error bars in their measured radio 
and X-ray luminosities thus implying a possible intrinsic scatter 
in the correlation. For this reason, as in Cassano et al.(2006), in 
our analysis we use the linear regression algorithm by Akritas & 
Bershady (1996) that indeed accounts for both intrinsic scatter 
and measured errors in both variables. The fits have been per- 
formed in the form : 

Log(Pi.4) -Y^A+b [Log(Lx) - X] (1) 

where P1.4 is in W/Hz, Lx is in erg/s, F=24.5, and X-A5 and 45.4 
in the case of the 0.1-2.4 keV luminosity (Fig. 1 Left) and bolo- 
metric X-ray luminosity (Fig. 1 Right), respectively. The best-fit 
normalizations and slopes of the correlations for giant radio ha- 
los, and the measured scatters across the correlation are given 
in Table 2. We report best fits obtained from BCES-bisector and 
orthogonal approaches and from their bootstraps (10000 boot- 
strap resamplings). The BCES bisector approach treats the vari- 
ables symmetrically and is recommended for scientific problems 
where the goal is to estimate relationships between the variables 
and for a comparison with theory (e.g., Isobe et al. 1990). 

Regardless of the nature of the scatter of the datapoint across 
the correlation, we point out that present data allow us to fairly 
constrain the slope of the correlations Q. 

3. Evolution of radio halos and connection with 
cluster mergers 

Correlations between radio and thermal properties in galaxy 
clusters may be explained by both hadronic and turbulent- 

' These slopes are consistent with those found in previous papers 
(e.g., Bacchi et al. 2003; Cassano et al. 2006). Kushnir et al.(2009) re- 
cently found significantly flatter slopes, however they do not fit the data 
with their errors measured in both variables but assuming an error in 
Pi 4 equal to the scatter of data points across the correlation. 
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Table 1. Radio and X-ray properties of clusters used in this paper: giant radio halos (first sector), mini-halos in cool core clusters 
(second sector), small scale halos (tirdth sector), clusters with no Mpc-scale radio emission (last sector). GMRT clusters are reported 
in bold. In Col.(l): Cluster name. Col.(2): Cluster redshift. Col.(3): X-ray luminosity in the energy range [0.1 - 2.4] keV in unit of 
10"*^ erg/s. Col. (4): Bolometric X-ray luminosity in the energy range [0.01 - 40] keV in unit of hj^ 10** erg/s. Col. (5): Radio 
power at 1.4 GHz in unit of 10^* Watt/Hz. Col.(6) References: 1 = Boehringer et al 2004, 2 = Ebeling et al 1998, 3 = Ebeling 
et al 1996, 4 = Tsuru et al 1996, 5 = Ebeling et al 2007, 6 = Ebeling et al 2000, 7 = Liang et al. 2000, 8 = Feretti et al. 2001, 9 
= Govoni et al. 2001, 10 = Bacchi et al 2003, 11 = Giovannini & Feretti 2000, 12 = Dallacasa et al 2009, 13 = Clarke & Ensslin 
2006, 14 = Brentjens 2008, 15 = Feretti 2002, 16 = Kim et al. 1990, 17 = Deiss et al. 1997, 18 = Govoni et al. 2005, 19 = Venturi 
et al 2007, 20 = Giacintucci et al 2009, 21 = van Weeren et al 2009, 22 = Bonafede et al 2009, 23 = Cassano et al 2008b, 24 = 
Giacintucci 2007 (integrating the diffuse radio emission between the two radio reUcs). (*) /'1.4 of A209 and RXJ1314 is estimated 
from that at 610 MHz by adopting a = 1.2. 



cluster's 


z 


Lx 


L'bol 


PlA 


Ref 


name 




[10^ erg/s ] 


[10^ erg/s ] 


[Vf^ Watt/Hz] 




1E50657-558 


0.2994 


23.03 ± 1.81 


87.10 ±6.82 


28.21 + 1.97 


1,7 


A2163 


0.2030 


23.17 ± 1.48 


86.50 ± 5.52 


18.44 ±0.24 


1,8 


A2744 


0.3080 


12.92 ± 2.41 


36.81 ±6.88 


17.16 ± 1.71 


1, 9 


A2219 


0.2280 


12.73 ± 1.37 


38.55 ±4.05 


12.23 ± 0.59 


2, 10 


CL0016+16 


0.5545 


18.83 ± 1.88 


55.59 ±5.56 


6.74 ± 0.67 


4, 11 


A1914 


0.1712 


10.71 ± 1.02 


34.51 ±3.28 


5.21 ± 0.24 


3, 10 


A665 


0.1816 


9.84 + 1.54 


27.61 ±4.32 


3.98 ± 0.39 


2, 11 


A520 


0.2010 


8.83 + 1.99 


23.77 ± 5.34 


3.91 ±0.39 


2, 10 


A521 


0.2475 


8.18 + 1.36 


19.63 ± 3.26 


1.16±0.11 


1, 12 


A2254 


0.1780 


4.32 ± 0.92 


11.35 ±2.43 


2.94 ± 0.29 


3, 10 


A2256 


0.0581 


3.81 ±0.17 


9.54 + 0.43 


0.68 ±0.12 


3, 13, 14 


A773 


0.2170 


8.10 ± 1.35 


22.70 ± 3.78 


1.73 ±0.17 


2,9 


A545 


0.1530 


5.66 ± 0.49 


12.50 ± 1.09 


1.48 ±0.06 


1, 10 


A2319 


0.0559 


7.40 + 0.40 


21.42 ± 1.16 


1.12±0.11 


3, 15 


A1300 


0.3071 


13.96 ± 2.05 


41.97 ±6.17 


6.09 ± 0.61 


1, 15 


A 1656 (Coma) 


0.0231 


3.77 ±0.10 


10.44 ± 0.28 


0.72+HI 

-0.04 


3, 16, 17 


A2255 


0.0808 


2.65 ±0.12 


6.61 ±0.30 


0.89 ± 0.05 


3, 18 


A754 


0.0535 


4.31 ±0.33 


12.94 ± 0.99 


1.08 ±0.06 


3, 10 


A209 


0.2060 


6.29 ± 0.65 


16.26 ± 1.69 


1.19 ±0.26 


1, 19* 


RXJ2003 


0.3171 


9.25 ± 1.53 


27.23 ± 4.95 


12.30 ± 0.71 


1,20 


MACS J0717 


0.5548 


24.6 ± 0.3 


84.18 ± 1.01 


50.0 ± 10.0 


5, 21, 22 


A2390 


0.228 


13.49 ±3.16 




9.77 ± 0.45 


2, 10 


Z7160 


0.258 


8.51 ±2.12 




2.19 ±0.26 


2, 23 


RXJ1314 


0.2439 


10.96 ± 1.81 




0.75 ± 0.15 


1, 19*, 24* 


A2697 


0.2320 


6.88 ±0.85 




<0.40 


1,25 


A141 


0.2300 


5.76 ± 0.90 




<0.36 


1,25 


A3088 


0.2537 


6.95 ± 1.20 




<0.42 


1,25 


RXJ1115.8 


0.3499 


13.58 ± 2.99 




<0.45 


1,25 


S780 


0.2357 


15.53 ±2.80 




<0.36 


1, 25 


RXJ1512.2 


0.3152 


10.19 ± 1.76 




<0.63 


1,25 


A2537 


0.2966 


10.17 ± 1.45 




<0.50 


1,25 


A2631 


0.2779 


7.57 ± 1.50 




<0.39 


1,25 


A2667 


0.2264 


13.65 ± 1.38 




<0.42 


1,25 


RXJ0027.6 


0.3649 


12.29 ± 3.88 




<0.68 


6, 25 


A611 


0.2880 


8.86 ± 2.53 




<0.40 


6, 25 


A781 


0.2984 


11.29 + 2.82 




<0.36 


2, 25 


Z2089 


0.2347 


6.79 + 1.76 




<0.27 


2, 25 


Z2701 


0.2140 


6.59 ± 1.15 




<0.42 


2, 25 


A1423 


0.2130 


6.19 ± 1.34 




<0.41 


2, 25 


Z5699 


0.3063 


8.96 ± 2.24 




<0.54 


6, 25 


Z5768 


0.2660 


7.47 ± 1.66 




<0.36 


6, 25 


Z7215 


0.2897 


7.34 ± 1.91 




<0.55 


6, 25 


RXJ1532.9 


0.3450 


16.49 ±4.50 




<0.62 


2, 25 


RXJ2228.6 


0.4177 


19.44 ± 5.55 




<0.91 


6, 25 



acceleration models with the observed slopes that can be repro- 
duced provided that also the magnetic field scales with cluster 
mass, temperature or luminosity (e.g. Dolag & Ensslin 2000; 
Miniati et al. 2001; Cassano et al. 2006; Dolag 2006). 



The new result from the high sensitivity data of the "GMRT 
Radio Halo Survey" is the bi-modal behaviour in Fig. 1, with 
radio-halo clusters and clusters without radio halos clearly sep- 
arated (Brunetti et al. 2007). It is not obvious to understand how 
two well separated classes of galaxy clusters can be generated 
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D 


C b 


0"Log(P,.4) 


O" LogiLx) 


p J 

^1.4~^0. 1-2.4 


RPFS Risertnr 


0.195 


0.060 


2.06 


0.20 


0.28 


0.14 




bootstrap 


v.Lyz 


n OAT 


Z.Uo 


u.zz 








BCES Orthogonal 


0.204 


0.062 


2.21 


0.23 








bootstrap 


0.202 


0.067 


2.23 


0.27 








BCES Bisector 


0.077 


0.057 


1.76 


0.16 


0.26 


0.15 




bootstrap 


0.075 


0.059 


1.78 


0.18 








BCES Orthogonal 


0.077 


0.059 


1.86 


0.19 








bootstrap 


0.075 


0.062 


1.88 


0.22 







Table 2. Fitting parameters for the correlations between f 1 4-L0.1-2.4 and Pi^-Lboi of giant radio halos (see text). 



in the context of the hierarchical process of large scale structure 
formation. Indeed clusters with similar thermal properties (mass, 
X-ray luminosity, . . . ) are expected to have a similar probability 
to host radio halos. In this case the observed differences in terms 
of cluster non-thermal properties should be understood by as- 
suming different evolutionary stages of these clusters. 
Clusters in the GMRT sample are selected with similar X-ray 
luminosity (simass) and redshift. The radio halo - merger con- 
nection (e.g. Venturi et al. 2008) may suggest that the difference 
between giant radio halo and "radio quiet" clusters is due to their 
dynamics. Thus regardless of the details of the mechanisms that 
generate radio halos, we shall assume the following evolutionary 
cycle : 

- /) galaxy clusters host giant radio halos for a period of time, 
in connection with cluster mergers, and populate the P\ 4-Lx 
coiTelation; 

- ;/) at later times, when clusters become dynamically relaxed, 
the Mpc-scale synchrotron emission is gradually suppressed 
and clusters populate the region of the upper limits. 

In this case, when restricting to clusters of the GMRT com- 
plete sample. Figure 1 provides a fair statistical sampling of 
the evolutionary flow of X-ray luminous clusters in the f 1 4-Lx 
plane at z=0.2-0.4. 

Radio-halo clusters, always dynamically disturbed systems, 
must be the "youngest" systems, where an ongoing merger, 
leading to their formation (or accretion of a sizable fraction of 
their mass), is still supplying energy to maintain the synchrotron 
emission. 

On the other hand, clusters with radio upper limits, typically 
more relaxed than radio halo clusters (Venturi et al. 2008), must 
have experienced the last merger at earlier epochs : after the last 
merger they already had sufficient time for suppression of the 
synchrotron emission and consequently they should be the "old- 
est" systems in the GMRT sample. 

Clusters in the "empty" region may be (a) "intermediate" sys- 
tems at late merging phases, where synchrotron emission is be- 
ing suppressed, or (b) the "very young" systems in the very 
early phases of a merging activity, where synchrotron emission 
is increasing. One cluster in the GMRT sample is found in the 
"empty" region, this is the merging cluster RXJ1314 that hosts 
a small-scale radio halo and 2 radio relics (Feretti et al. 2005; 
Venturi et al. 2007) and consequently it likely belongs to the lat- 
ter class (b) of galaxy clusters (although we cannot exclude that 
it is an "intermediate" system). 

4. Constraining the evolution of radio halos 

In the case that we admit that radio halos are transient phe- 
nomena connected with cluster-merging phases, the "emptiness" 



of the region between radio halos and "radio quiet" clusters in 
the P\ 4-Lx diagram can be used to constrain the time-scale of 
the evolution (suppression and amplification) of the synchrotron 
emission in these clusters (Brunetti et al. 2007). Indeed the sig- 
nificant lack of clusters in this region suggests that this time- 
scale is much shorter than both the "life-time" of clusters in the 
sample and the period of time clusters spent in the radio halo 
stage. 

By restricting our analysis to the 19 clusters in the GMRT sam- 
ple with Lx ^ 8.5 X 10"*"* erg s"*, in which case the radio power of 
giant radio halos is 1 order of magnitude larger than radio upper 
limits, we find 5 giant radio halos (and 2 mini-halos in cool-core 
clusters) on the coiTelation, 1 1 clusters in the region of the up- 
per limits and only RXJ1314 in the "empty" region. Thus, the 
time interval that clusters may spend crossing this "empty" re- 
gion is Si fTgc, where / is the fraction of clusters in this region, 
1/19, and Tgc is the period of time elapsed since the last merger 
in the case of the "oldest" clusters in our population. Since the 
GMRT sample is constituted by massive, M > IO'^Mq, galaxy 
clusters at z=0.2-0.4, Tgc is essentially the time-scale between 
the epoch of formation of these clusters (z»0.6-0.7, e.g. Giocoli 
et al. 2007) and the most recent epoch of observation (z=0.2, 
w 11.2 Gyr), t^c ~ 3.5 Gyr. Consequently, the "life-time" of 
radio halos is Trh = 7/19Tgc w 1.3 Gyr, the time interval that 
clusters may spend in the "empty" region is = 1/19ti,c ~ 180 
Myr, and the coiTesponding time-scale for suppression of the 
cluster-scale synchrotron emission from the level of radio halos 
to that of "radio quiet" clusters, r, is roughly half of this period, 
T w 90 Myrs, considering that clusters may cross the "empty" 
region two times (during the amplification and suppression of 
synchrotron emission. Sect. 3). We stress that this conclusion 
holds even if the two mini-halos, Abell 2390 and Z7160, are ex- 
cluded from our analysifl in this case t,.^ w 1 Gyr, / - 1/17, 
and T K 200 Myrs. 

Obviously the poor statistics allows fairly large uncertain- 
ties on the above numbers. Consequently, because in Sect.5 we 
will discuss the implications of these constraints for the origin 
of radio halos, it is important to understand whether, due to the 
poor statistics, a significantly larger value of t is still consis- 
tent with the distribution of galaxy clusters in Figure 1 . We use 
Montecarlo procedures : we assume that clusters spend a frac- 
tion of time fi„„, flip and / in the region of the upper-limits, 
on the coiTelation and in the "empty" region, respectively. Then 
we perform 10^ random extractions of N synthetic clusters with 
probability /,„„., fup and / (with /,„„. + fup + f = 1), and de- 
rive the fraction of trials where 1 cluster falls in the "empty" 



^ Mini-halos in cooling-flow clusters may have a different origin with 
respect to giant radio halos, possibly connected with the presence of the 
cooling flow (Gitti et al. 2002; Mazzotta & Giacintucci 2008; Cassano 
et al. 2008b). 
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Fig. 2. Left panel: fraction of trials from Montecarlo simulations that match observations (i.e. 1 cluster in the "empty" region) as 
a function of / (see text), assuming N= 17 (solid line), 19 (dotted line) and 20 (dashed line) (Table 3). Right panel: distributions 
of trials from Montecarlo simulations as a function of the number of clusters in the empty region, N^mp, assuming /=0.125 (solid 
line and filled circles), 0.17 (dotted line and empty circles), 0.3 (dashed line and empty squares), 0.4 (long-dashed line and empty 
triangles), r - 3.5/ and = 3.5//2 Gyrs if one assumes that clusters cross the "empty" region one or two times, respectively. The 
vertical dashed line marks 1 cluster in the "empty" region. 



region. This is shown in Figure 2 (Left) as a function of / : 
as expected the distribution peaks at / ^ 0.06, corresponding 
to T ~ 1/2/Tgc ~ 100 Myr, in which case about 38% of trials 
match observations. In the less constrained case, N-ll, where 
the 2 mini-halos in the GMRT sample are not considered, the 
fraction of trials that match observations falls to only 1 % and 
0.3 % for /=0.33 (t = 0.57) and /=0.39 (0.68 Gyrs), respec- 
tively (see Table 3, where we also report the results obtained by 
considering different sub-samples). 

This analysis allows us to conclude that values of t signif- 
icantly larger than a few tenths of Gyr are very unlikely. This 
is also highlithed by Figure 2 (Right) that shows the distribu- 
tion of trials (in the less constrained case, N-\l) as a function 
of the number of galaxy clusters found in the "empty" region 
for different values of t. Larger values of t imply an increasing 
number of clusters expected in the "empty" region, and t > 0.6 
Gyr can be excluded at > 99% confidence level. We note that our 
conclusion is inconsistent with much larger values of the transi- 
tion time-scale, r ~ few Gyrs, as recently claimed by Kushnir 
et al. (2009) that however estimate the transition time-scale as 
T ~ l/3Ti;c, 1/3 being the fraction of clusters with radio halos 
and Tgc taken = 5 Gyfl On the other hand, our statistical anal- 
ysis provides more quantitative support to previous conclusions 
(Brunetti et al. 2007; Brunetti 2008). 

For completeness, we also consider the complementary sce- 
nario where clusters cross the "empty" region only one time, due 
to the suppression of their synchrotron emission Q. In this case 
we interpret RXJ1314 as an "intermediate" system in early -post 
merging phase (Sect. 3) that provides the most conservative ap- 



^ This approach indeed would give the life-time of radio halos, not 
the transition time-scale, and it is indeed consistent with our estimate of 

Trh- 

^ This might happen if the increase of the cluster X-ray luminosity 
during mergers takes longer times than that of the synchrotron luminos- 
ity, and consequently merging clusters may approach the range of X-ray 
luminosities of the GMRT sample "along" the correlation. 



proach to constrain r (Table 3). Still also in this conservative 
approach we conclude that present data strongly favour values 
of T substantially smaller than 1 Gyr (r > 1 Gyr is excluded at 
> 99% in our reference case, N - 19, Table 3). 

5. Implications for the origin of giant radio halos 

5.1. Hadronic models 

Theoretically relativistic protons are expected to be the dominant 
non-thermal particle component in galaxy clusters since they 
have very long life-times and remain confined within clusters for 
an Hubble time (e.g. Blasi et al. 2007 and ref. therein). Proton- 
proton (p-p) collisions provide a continuous source of secondary 
products in the ICM, and secondary electrons in turns generate 
diffuse synchrotron emission. 

Radio halos are found in merging clusters and the passage of 
merger shocks through the ICM may increase the energy density 
of protons (e.g. due to acceleration of these protons at merger 
shocks) enhancing the rate of production of secondary electrons 
and the resulting cluster-scale synchrotron emission. However, 
since protons have very long life-times the production rate of 
secondary electrons in the ICM would remain unchanged with 
cosmic time and the mechanism itself does not allow a suppres- 
sion of the synchrotron emission when clusters become more 
dynamically relaxed. 

Clusters of equal masses may experience different formation his- 
tories, yet the global budget of gravitational energy dissipated at 
merging and accretion shocks is expected to be similar yielding 
to fairly small diff'erences in terms of energy content of rela- 
tivistic protons (Jubelgas et al. 2008). Consequently, assuming 
that the ICM is magnetised at fiG level, radio halos generated 
by secondary emission are expected long-living and common; 
also, some trend between their radio power and the X-ray lumi- 
nosity or temperature of the hosting clusters is expected (Dolag 
& Ensslin 2000; Miniati et al. 2001; Dolag 2006; Pfrommer et 
al. 2008). Due the huge uncertainties in the physics of shock 
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Table 3. Synchrotron dissipation time (in unit of Gyr) in galaxy clusters that allows to match observations (1 cluster in the 'empty" 
region) in a fraction of Montecarlo trials = 10% (col.5), 5% (col.6), 1% (col.7) and 0.3 % (col.8). Montecarlo simulations are carried 
out for three configurations: GMRT clusters with Lo.i-2.4 > 8.4 ■ lO^'^erg s"' excluding the two mini-halos (first line), GMRT clusters 
with L().i-2.4 > 8.4- 10'*'*erg s"' (second line), GMRT clusters withL() i_2.4 > 8 ■ 10'*'*erg s"' excluding the two mini-halos (tirdth line); 
total number of observed clusters and their distribution are given in col. 1-4. We perform 10^ random extractions of TV synthetic 
clusters assuming a grid of extraction probability with / in the range 0-1 (see text). Upper part gives the case where clusters cross 
the "empty" region two times, lower part gives the case where clusters cross this region only one time. 



acceleration one may easily believe that large variations of the 
content of relativistic protons are possible among clusters with 
similar mass (temperature ..), however there is no reason to ex- 
pect a bi-modality in the cluster synchrotron emission. 

Consequently, to explain the separation between radio-halo 
and "radio quiet" clusters and the merger-halo connection the 
magnetic field should play a major role and it must be admitted 
that merging clusters, hosting radio halos, have larger magnetic 
fields and that this excess in magnetic field is dissipated when 
clusters become "radio quiet" and dynamically relaxed (Brunetti 
et al. 2007, 2008; Kushnir et al. 2009). 

Synchrotron emission in hadronic models scales as (e.g. Dolag 
& Ensslin 2000): 



(2) 



where Z?„„/, = 3.2(1 +z)^juG is the equivalent field due to inverse 
Compton scattering of Cosmic Microwave Background photons 
and a ~ 1.3 is the synchrotron spectral index of radio halos 
(e.g. Ferrari et al. 2008). Thus to explain a suppression >10 in 
terms of synchrotron emission (Figure 1) the ratio between the 
magnetic fields in radio halos, B + 6B, and that in "radio quiet" 
clusters, B, must be : 
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Fig. 3. Lower limit to the ratio between the energy density of 
the magnetic field in radio halos and in "radio quiet" clusters 
as a function of the magnetic field strength in radio halos. The 
vertical dashed line marks the value of the equivalent field of the 
Cosmic Microwave Background photons assuming z=0.25. 
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The ratio between the magnetic field energy densities in the two 
cluster populations, (jJrhl(Orq <^ (1 +6B/B)^, fromEq. 3 is shown 
in Figure 3 for z ~ 0.25, typical of GMRT clusters. In the case 
B + 6B « Bcmb, hadronic models must admit that the energy 
density of the magnetic field in "radio quiet" clusters is >10 
times smaller than that in radio halos, and even larger ratios must 
be admitted in the case B + 6B >> B^j^ij. 

Theoretically we might admit that the magnetic field is am- 
plified in the ICM by turbulence generated in cluster mergers 
(Dolag et al. 1999, 2005; Ryu et al. 2008), and later dissipated 
since turbulent magnetic fields can decay. 
On the other hand, to our knowledge, studies of Faraday 
Rotation in galaxy clusters do not find any statistical difference, 
in terms of energy density of the large scale (10-100 kpc coher- 
ent scales) magnetic field, between clusters hosting radio halos 



and clusters without Mpc-scale radio emission (e.g. Carilli & 
Taylor 20020). 

Most important, dissipation of this magnetic field is expected 
to take long time. Even if we simply consider the case where 
the field is dissipated through the decay of cluster-MHD turbu- 
lence, the energy density of the rms field decreases only (about) 
linearly with the eddy turnover time-scale. This requires sev- 
eral eddy turnover times, « a few Gyr, to gradually dissipate 
the bulk (i.e. 80-90 %) of the energy density of the field in the 
ICM (Subramanian et al. 2006), that indeed also explains why 
few yuG-fields are common in galaxy clusters. Consequently the 
dissipation time-scale of the magnetic field is inconsistent with 
(larger than) that of the suppression of the cluster-scale syn- 
chrotron emission inferred from the statistical analysis in pre- 
vious Section. 



^ We would also point out that in some cases the magnetic field in 
"radio quiet" clusters, e.g. A119, is larger than that of radio halo clus- 
ters, A2255, with similar X-ray luminosity (Murgia et al. 2004, Govoni 
et al. 2005). 
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Fig. 4. The intensity of the small-scale magnetic field in radio 
halos that must be dissipated to match observations is reported 
as a function of the strength of the large-scale magnetic field in 
"radio quiet" (and radio halo) clusters. 



This conclusion is based on the scenario, following Subramanian 
et al. (2006), that the magnetic field in the ICM is amplified by 
cluster turbulence generated on large scale, ~ 150-300 kpc, in 
which case the thickness of magnetic filaments is expected to 
be ~20-40 kpc. This is supported by Faraday Rotation measure- 
ments of extended sources in clusters that allow to observe order- 
ing scales of the magnetic field ~ 10-40 kpc (Clarke et al. 2001; 
Guidetti et al. 2008) and that indicate, at least in some cases, 
that the power spectrum of the magnetic field extends to very 
large scales, 100-500 kpc (e.g. Murgia et al. 2004). On the other 
hand our understanding of the origin and of the properties of 
magnetic field in galaxy clusters is still poor and leaves space 
to large uncertainties. Indeed, faster dissipation of the magnetic 
field in galaxy clusters, r « several 100 Myr, may happen in the 
case that the magnetic field in excess, SB, in clusters hosting ra- 
dio halos is associated with a field component on smaller scales. 
The value of the small scale field, 6B, necessary to account for 
the difference between the synchrotron emission in radio halo 
and "radio quiet" clusters can be obtained from Eq. 3 and is re- 
ported in Figure 4 as a function of the large scale magnetic field, 
B. Figure 4 clearly highlights the drawbacks of this hypothesis 
: first of all the small scale field must be energetically dominant 
with respect to that on larger scales (see also Figure 3), in addi- 
tion if the large scale field is > 1.5 - 2jiG level (consistent with 
present RM studies) then 6B would be extremely large, > 10 /iG. 
Since 6B must be dissipated in a few tenths of Gyr, we note that 
this would imply a magnetic-energy dissipation rate in a Mpc^ 
region > 10^^(t /0.5Gyr)~^ erg/s, e.g. larger than the bolometric 
X-ray emission of clusters themselves. 

5.2. Turbulent acceleration of particles 

MHD turbulence generated during cluster mergers may accel- 
erate relativistic particles (e.g. Brunetti et al. 2008). Even with- 
out considering the dissipation (or amphfication) of the mag- 
netic field in the ICM, the finite dissipation time-scale of turbu- 
lence implies that giant radio halos should be found in merging- 
clusters where turbulence is still generated and must be ex- 




8.5 9 9.5 

Log(y[Hz]) 

Fig. 5. Example of the synchrotron emitted power (arbitrary 
units) as a function of frequency for different energy densities of 
turbulence (magnetosonic waves) : 30 % (sohd fine), 25 % (dot- 
ted Une), 20 % (dashed line), 15 % (long-dashed line) of the ther- 
mal energy density. In the calculations we adopt a homogeneus 
model with magnetic field strength - 3 //G, number density of 
the thermal plasma - lO^^cm"^ and temperature T-1 ■ lO^K. 
Frequency ranges of interest for GMRT, VLA and LOFAR are 
also marked. 



tremely rare in more relaxed clusters (e.g., Cassano & Brunetti 
2005). 

As soon as large scale turbulence in the ICM reaches smaller, 
resonant, scales (via cascading or induced plasma instabilities, 
e.g. Brunetti et al. 2004, Lazarian & Beresnyak 2006, Brunetti 
& Lazarian 2007), particles are accelerated and generate syn- 
chrotron emission. In the case of radio halos emitting at GHz 
frequencies the acceleration process should be relatively effi- 
cient and particles get accelerated to the energies necessary to 
produce synchrotron GHz-emission within a time-scale smaller 
than a couple of cooling times of these electrons, that is ~ 100 
Myrs. Although the large uncertainties in the way large scale 
turbulence is generated in the ICM during cluster mergers, it is 
likely that the process persists for a few crossing times of the 
cluster-core regions, that is fairly consistent with a radio halo 
life-time Trh ~ 1 Gyr as derived in Section 4. 
Most important, the cooling time of the emitting electrons is 
smaller than (or comparable to) the cascading time- scale of the 
large-scale turbulence implying that the evolution of the syn- 
chrotron power depends very much on the level of MHD tur- 
bulence in the ICM (e.g., Cassano & Brunetti 2005; Brunetti & 
Lazarian 2007). Consequently, if we simply assume that the in- 
jection of MHD turbulence is suppressed "instantaneously" at a 
given time (eg. at late merging-phase), then also the synchrotron 
emission at higher radio frequencies is suppressed, falling below 
the detection hmit of radio observations, as soon as the energy 
density of turbulence starts decreasing. This is shown in Figure 5 
where the synchrotron spectrum from turbulent accelerated elec- 
trons is reported for different energy density of the MHD turbu- 
lence. A reduction of the turbulent energy density of a factor 2 
happens within about 1 eddy turnover time of the large scale tur- 
bulence, that is a few times 100 Myrs, and this is sufficient to 
suppress the synchrotron emission at higher, GHz, frequencies 
by about 1 order of magntitude. 
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Consequently cluster bi-modality in this scanario may be ex- 
pected because the transition between radio halos and "radio 
quiet" clusters in the P\a-Lx diagram is expected to be fairly 
fast (Brunetti et al. 2007, 08) provided that the acceleration pro- 
cess we are looking in these sources is not very efficient, being 
just enough to generate radio halos emitting at a few GHz fre- 
quencies. On the other way round, we might say that the ob- 
served cluster bi-modality constrains the efficiency of the parti- 
cle acceleration process in radio halos. Interestingly, a relatively 
inefficient electron acceleration process in radio halos is in line 
with the steep spectrum observed in these sources and, most im- 
portant, with the presence of a spectral steepening at higher fre- 
quencies discovered in a few halos (e.g., Thierbach et al. 2003, 
Brunetti et al. 2008, Dallacasa et al. 2009). 
Also, Figure 5 suggests that under these conditions the sup- 
pression of synchrotron emission that follows the dissipation of 
MHD turbulence is more efficient at higher frequencies and thus 
cluster bi-modality is expected to be less pronounced in consid- 
ering the synchrotron emission of galaxy clusters at lower fre- 
quencies. This is a clear expectation of the scenario that can be 
tested with future observations of samples of galaxy clusters at 
100-200 MHz that may be carried out with LOFAR in a couple 
of years. 

6. Conclusion 

The "GMRT Radio Halo Survey" allows to study the statistics 
of radio halos in a complete sample of X-ray luminous galaxy 
clusters (Venturi et al. 2008). The high sensitivity of the ra- 
dio observations at the GMRT allows to unveil a cluster radio 
bi-modality with "radio quiet" clusters well separated from the 
region of the Pi^-Lx correlation defined by giant radio halos 
(Brunetti et al. 2007 and Figure 1). 

In the framework of the hierarchical model galaxy clusters 
are expected to evolve in the Pia-Lx plane, in which case the 
distribution of GMRT clusters in Figure 1 results from a sta- 
tistical sampling of this evolution. The connection between ra- 
dio halos and cluster mergers suggests that the Mpc-scale syn- 
chrotron emission in galaxy clusters is amplified during these 
mergers and then suppressed when clusters become more dy- 
namically relaxed. The separation between radio halo and "radio 
quiet" clusters in Figure 1, and the rarity of galaxy clusters with 
intermediate radio power suggests that the processes of ampUfi- 
cation and suppression of the synchrotron emission takes place 
in a relatively short time-scale. 

The time-scale of the evolution from radio halos to "radio 
quiet" clusters (and vice versa) provides a novel tool to con- 
strain models proposed for the origin of radio halos, namely the 
re-acceleration and hadronic model. In the former case the ac- 
celeration and cooling of relativistic electrons drive the level of 
the Mpc-scale synchrotron emission from clusters, while in the 
latter case the transition between radio halo and "radio quiet" 
clusters must be due to the amplification and dissipation of the 
magnetic field in the ICM. 

We carried out statistical analysis of the cluster radio bi-modality 
in Figure 1 and, although the still poor statistics, show that the 
suppression of the cluster-scale synchrotron emission must hap- 
pen in a fairly short time-scale, a few 100 Myrs, whereas longer 
time-scales, Gyr, are not consistent with present data. This short 
transition time-scale can be potentially reconciled with the hy- 
pothesis that the emitting electrons are accelerated by cluster- 
scale turbulence, in which case the synchrotron radiation emit- 
ted at GHz frequencies may rapidly decrease as a consequence 
of the dissipation of a sizeable fraction of that turbulence. In this 



case, however, we also claim that a less pronounced bi-modality 
is expected in the case of cluster samples observed at lower ra- 
dio frequencies, that may be tested by future LOFAR and LWA 

observations. 

On the other hand, it is more difficult to reconcile a short tran- 
sition time-scale in ffie case that the unique source of emitting 
electrons is provided by p-p colhsions (hadroiuc models). In this 
case ffie dissipation of the cluster magnetic field that suppresses 
the synchrotron emission would take longer periods of time. In 
principle this difficulty could be considerably alleviated in ffie 
case that ffie energy density of ffie magnetic field in radio halos is 
dominated by that of small scale field. However, we would come 
into ffie untenable scenario in which a very strong, transient mag- 
netic field (small scale) component, is present in the ICM. Future 
studies of source-depolarization in cluster radio sources will also 
help in constraining the level of the small-scale field component 
in ffie ICM. 

We stress that our constraints come from ffie conservative (and 
simplified) assumption that the injection of turbulence (as well 
as the amplification of the magnetic field) switches ofi" at the 
same time across the radio halo, Mpc^, region. In reality, de- 
pending on the way turbulence and large scale magnetic fields 
are generated in ffie ICM, the suppression and amplification of 
the synchrotron emission could start at different times in diflFer- 
ent parts of this volume. The most important consequence of 
that is an expected scatter in ffie correlation raffier ffian in the 
way clusters become "radio quiet", since clusters are expected to 
start moving across ffie transition region only when synchrotron 
is suppressed across a substantial fraction of ffie radio halo's vol- 
ume. Yet, overall this goes into the direction to strengffien our 
conclusion that an efficient process to suppress the cluster-scale 
synchrotron emission in galaxy clusters is necessary to explain 
observations. 

Giant radio halos prove complex physical processes where a 
fraction of the gravitational energy dissipated during cluster- 
mergers is channelled into ffie acceleration of relativistic parti- 
cles. The correlation traced by halos in Figure 1 and its intrinsic 
scatter, together wiffi the disttibution of clusters in ffie Pi^-Lx 
plane, provide novel tools to hopefully constrain the complex 
physics of turbulence and magnetic fields in the ICM and their 
interplay wiffi ffie process of cluster formation. The deep sur- 
veys at low frequencies with LOFAR and LWA will be crucial to 
overcome present uncertainties due to ffie still poor statistics al- 
lowing a major step forward in understanding ffie origin and evo- 
lution of the cluster-scale synchrotron emission. Remarkably, as 
discussed in this paper, the scatter of the correlation and ffie dis- 
tribution of clusters in ffie Pi^-Lx plane are expected to depend 
on ffie observing radio frequency (Sect. 5.2, Figure 5), and con- 
sequently deep complementary follow ups at intermediate and 
higher frequencies (GMRT, eVLA, SKA) will also be crucial. 
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